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Abstract. - We have investigated the proximity effect between a superconductor (Nb) and 
a 'Heavy Fermion' system (CeCue) by measuring criticai temperatures T c and paraliei criticai 
fiefds _ffj.' 2 (T) of Nb fiims with varying thickness deposited on 75 nm thick films of CeCug, and 
comparing the results with the behavior of similar films deposited on the normal metal Cu. For 
Nb on CeCu6 we find a strong decrease of T c with decreasing Nb thickness and a finite critical 
thickness of the order of 10 nm. Also, dimensional crossovers in i?| 2 (T) are completely absent, 
in strong contrast with Nb/Cu. Analysis of the data by a proximity effect model based on the 
Takahashi-Tachiki theory shows that the data can be explained by taking into account both 
the high effective mass (or low electronic diffusion constant), and the large density of states at 
the Fermi energy which characterize the Heavy Fermion metal. 



Introduction. - When a superconductor (S) is brought into contact with a non-superconducting 
conductor (N) , superconductivity leaks into that material by the proximity effect pQ . For the 
thermodynamic properties of the S /N bilayer such as the critical temperature T c or the upper 
critical field H C 2 , this can be modelled as a spatial variation of the superconducting pair den- 
sity. The spatial variation mainly depends on the electron diffusion constants for the S- and 
the N-metal, on the transparency of the interface for electrons and Cooper pairs, and on the 
pair breaking mechanisms on the N-side of the interface. Until now, two classes of N-materials 
have been investigated. The first is formed by simple metals such as Cu or Au, where the 
physics is mostly understood. Pairs are broken at finite temperatures by thermal fluctuations, 
leading to dephasing of the constituents of the induced Cooper pair. This is translated into 
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a temperature-dependent characteristic length over which superconducting correlations pen- 
etrate, the 'normal metal coherence length' £jy Since £jv = \J (frD^)/ ftirksT) {Dm is the 
diffusion constant of the N-metal, T is the temperature, the other symbols have their usual 
meaning), this length can become large at low temperatures. In the absence of other pair 
breakers, T c of a bilayer with finite N-layer thickness <In will be finite. The other class of 
materials is formed by ferromagnets (F) such as Fe; pair breaking is due to the exchange in- 
teraction E ex which acts on the spins of the Cooper pair. For strong magnets (E ex » ksT c ) 
it results in a temperature-independent £p = \J (TiDp) / (2irE ex ) with a characteristic value 
of only a few nanometer. In this case, superconductivity can already be quenched at a finite 
value for ds, called the critical thickness d c g 

Little attention has yet been paid to N-layers where the electronic ground state is dominated 
by many body correlations, such as in a Heavy Fermion (HF) metal. Basically, HF metals 
consist of a lattice of atoms with localised (f-)electrons, where the magnetism is quenched 
by a coherent Kondo effect. This leads to a strong peak in the DOS near the Fermi en- 
ergy with small energy width (for relevant reviews, see The consequences for the low 
temperature physics can be phenomenologically described in terms of Landau Fermi liquid 
theory by a mass renormalization of the charge carriers. The ensuing large effective mass m* 
can be directly observed in e.g. the specific heat of the system, c v , where the linear term 
7 = Cv/T oc m* /m e (with m e the bare electronic mass) can be up to two orders of magnitude 
larger than in normal metals. Since the magnetic moments are often not completely quenched, 
HF systems can also show magnetic order at low temperatures. For the proximity effect in an 
S/HF system, several scenarios are possible. In the spirit of the Fermi liquid approach, the 
HF metal can be considered a normal metal with a large mass, or equivalently a low Fermi 
velocity and therefore a small diffusion constant. Also, additional pair breaking mechanisms 
may be present due to the strong electron-electron interactions or the residual moments. On 
the other hand, the interface transparency may be decreased due to the mismatch in Fermi 
velocities, which would shield the superconductor from the HF metal and counter the other 
two effects. 

In order to investigate these questions, we have performed a comparative proximity effect 
study of the thin film systems Nb/CeCu6 and Nb/Cu. CeCu6 is a well-known HF system, 
with an extremely large value of 7 « 1.6 J/(mole K 2 ) j^] and no magnetic order down to the 
mK- regime, making it a good model system for this study. Nb/Cu is a much studied prox- 
imity system which shows one feature of particular interest, namely a Dimensional Crossover 
(DCO) from three-dimensional (3D) to two-dimensional (2D) behavior in the temperature 
dependence of the critical field parallel to the layers h\ 2 HE]- The DCO occurs when the 
N-layer thickness becomes of the order of and can be modelled quite accurately. In par- 
ticular, recent calculations by Ciuhu and Lodder based on the Takahashi-Tachiki theory for 
metallic multilayers and including a finite interface resistance Rb, showed that quantita- 
tive agreement between theory and experiment can be found in the case of Nb/Cu for very 

reasonable values of the different parameters Here we present data on T c {ds) and H\ 2 {T) 
for bilayers Nb/Cu and Nb/CeCu6 and compare them with similar calculations. We show that 
the HF system can be treated as a normal metal, with both the low diffusion constant (Fermi 
velocity) and the high DOS as necessary ingredients to explain the observed behavior. 

experimental. - Sets of CeCug films and bilayers of swfe/CeCug/Nb {sub denotes the sub- 
strate) were fabricated by DC-magnetron sputtering in an ultra high vacuum (UHV) system 
with a background pressure of the order of 5 x 10~ 10 mbar, and an Ar sputtering pressure of 
2.5 x 10~ 3 mbar. Crystalline CeCu6 was grown as reported before at a temperature of 
350 °C, using Si-substrates with amorphous SisN4 buffer layers to prevent Cu diffusion at those 
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Fig. 1 - Materials characteristics of the CeCu6 films, (a) RBS spectrum (nr. of counts versus 
backscatter energy Ess) taken with 4 He-ions of 2 MeV on a sample Si/Si3N4/(75 nm CeCue)/(15 
nm Nb). The different elements are indicated. The thin smooth line is a fit to the measured curve, 
(b) Specific resistance as function of temperature for single CeCug films of 75 nm (o) and 100 nm 
(A) . 

temperatures. The Nb was deposited on top of the CeCu6 after cooling the substrate holder 
with cold nitrogen gas to close to room temperature. Composition, thickness and crystallinity 
of the films were determined by Rutherford BackScattering (RBS) measurements together 
with X-ray diffraction measurements at low and high angles. The RBS measurements show 
good agreement with the expected stoichiometry for CeCug and no diffusion is found either of 
Ce or Cu into the substrate or of Nb into the CeCug. Fig.^ shows part of the RBS spectrum 
for CeCug(75 nm) / Nb(15 nm) on a Si/SiaN4 substrate and a fit of the data without taking 
any diffusion into account. Bilayers and trilayers of sit&/Cu/Nb and Cu/Nb/Cu were grown 
in a different UHV system with similar background pressure and sputtering conditions. In 
order to compare results, d,Nb in the trilayers was taken two times d,Nb in the bilayers, which 
yields equal conditions for the superconducting order parameter in the middle of the film (for 
the trilayer) and at the vacuum interface (for the bilayer). 

Resistance measurements were performed in standard 4-point geometry on lithographically 
patterned samples with bridge widths of 200 fim and a distance of 1.2 mm between the volt- 
age contacts. The electrical resistivity pceCue as function of temperature for two single films 
is shown in Fig. ^> and behaves as reported before |llj , with a clear maximum in i?(T) at 
T mal ss 5 K, similar to what is found for bulk material and a residual resistivity of 
the order of 80 /jficm. We decided to use 75 nm thick CeCug layers, which with T max =A K 
suggest only little deviation from the bulk properties. 

Results and discussion. - The dependence of T c on djv& for the bilayer set s«6/CeCu6(75 
nm)/Nb(djvb) is shown in Fig.[3i and compared to that of single Nb films, prepared in both 
UHV systems. For the Nb films a slight suppression is witnessed, usually found in the case 
of Nb, and caused by a mixture of different effects such as oxidation through the unprotected 
top layer, a small proximity effect with the substrate, and lifetime broadening due to growth 
disorder ^3]. For the S/HF bilayers, the suppression is much stronger, with a critical thickness 
d cr for onset of superconductivity reached around 12 nm. This would be equivalent to 24 nm in 
a trilayer configuration, and is of a similar magnitude as found in superconductor/ferromagnet 
systems such as V/Fe and Nb/Fe |3lll3) . The result shows immediately that the interface does 
allow particle exchange, so the supposed huge Fermi velocity mismatch in the system of order 
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Fig. 2 - (a) T c as function of Nb thickness d^b for bilayers sub/ 'CeCu6 /Nb (»)and single Nb films 
(+). The drawn line is a fit using the proximity effect model. The dotted line is a model calculation 
with NceCud = Ncu- The inset shows the behavior of -Djvt as found from the fit. (b) Parallel critical 
field H\ 2 as function of T for a bilayer sit&/Cu/Nb (A;djv6 = 15 nm) and a trilayer Cu/Nb/Cu (o; 
djv6 = 30 nm). Drawn lines are model fits, (c) h} 2 as function of T for bilayers sufe/CeCue/Nb with 
thicknesses djvt = 13.5 nm (lowest T c ), 15 nm, 17 nm, 25 nm, 30 nm, 50 nm. Drawn lines are model 
fits, (d) Model calculation of the inverse slope of Hj 2 at T versus the effective electron mass m*. 



10 3 does not lead to a highly reflective interface. It also shows that the coherence length in 
the HF-material ^hf must be small, leading to the strong suppression. 

Another indication of behavior deviating from simple S/N systems comes from the parallel 
critical field tljJ 2 (T). To demonstrate the difference, Fig. |3* shows h| 2 (T) for a bilayer of 
sub/ Cu(75 nm) / Nb(15 nm) and a trilayer of Cu(75 nm) / Nb(30 nm) / Cu(75 nm), to be 
compared to the data of the S/HF bilayers given in Fig. |2t- Decreasing T near T c , hJJ 2 (T) is 
linear for the Nb/Cu samples, followed by a kink and wl — T/T c -like 2D behavior. The kink 
signals the well-known DCO |7||^j. usually observed in multilayers, but also present in tri- and 
bilayers. In strong contrast, the S/HF bilayers do not show a DCO but only 2D-behavior for 
all dpfb- Qualitatively, this again indicates a small value for : the superconducting order 
parameter does not penetrate sufficiently far into the HF metal to yield a coupled system. 
To make more quantitative statements, we analyzed the data by means of model calculations 
based on the Takahashi-Tachiki formalism. Details are given in ref. [TJ3]) but we reiterate 
the main elements in order to introduce the different parameters. The formalism solves the 
equation for the pair potential A(r) with a space-dependent coupling constant V(r) for small 
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metal 


d,N 


D N 




Rb 


Cu (bi) 


15 


240 


1.4 


177 


Cu (tri) 


30 


220 


2.9 


252 


CeCu 6 


13.5 


0.15 


1.2 


324 




15 


0.19 


1.4 


338 




17 


0.11 


1.2 


350 




25 


0.11 


1.75 


252 




30 


0.11 


2. 


252 




50 


0.11 


2.2 


324 



Table I - Values for the normal metal thickness d,N (in nm), the fitted values for the diffusion 
constants Dcu, DceCue, Dm (in 10~ 4 m 2 / s) and the interface resistance Rb (in 10~ 8 /iQcm 2 ) for 
the different samples discussed in the text 



A close to T c : 

A(r) = V(r)k B Tj2 I d 3 r'Q u (r, r')A(r') (1) 

in which Q u is an integration kernel still to be determined, V(r) is the BCS coupling constant 
and the other symbols have their usual meaning; the summation runs over the Matsubara 
frequencies. In the dirty limit, can be shown to obey : 

[2| W |~^(r)(V--^A(r)) 2 ]Q w (r,r') = 2kN{t)8{t - r') (2) 
Tic 

with N(r) the DOS at the Fermi energy and D(r) the diffusion constant. These equations are 
complemented with boundary conditions for A at the interface which parametrize a possible 
barrier encountered by the electrons through a boundary resistivity Rb ■ In the calculations we 
use fixed values for the density of states ratios Ncu/N^b = 0.16 and NceCue/^Nb — 320. The 
former value was also used in ref . ^U] an d derives from ref . ^3] ; the latter value is constructed 
from the former by the ratio JCeCue/lCu = 2000. The value of T Cj Nb was fixed at 9.2 K, 
except for the Nb/Cu bilayer with d,Nb = 15 nm where it was 8.2 K. This reflects the fact that 
T c for thin Nb-layers starts to decrease, as explained above. Fitted were the different diffusion 
constants Dn and boundary resistances Rb- The results of the fits for the critical field data 
are shown in Fig. [21 as solid lines. The parameter values are given in Table |I| For the two 
Nb/Cu samples, the fitted values are in very reasonable agreement with the numbers found 
by fitting the data of Chun et al [SlllOj. The values for D^f, are roughly equal, the values for 
Dcu are slightly higher, which probably reflects the difference in preparation conditions, and 
also the values for Rb are very similar For the Nb/CeCu6 samples it can be noted that 
Djvb increases slowly and more or less linearly with increasing d^b] using this linear variation, 
which is plotted in the inset of Fig. 03, and values of DceCu6 = 0.1 xl0 _4 m 2 s~ 1 , Rb = 324 
x 10 -8 /ificTO 2 we calculated the behavior of T c (ds) as a consistency check. The agreement, 
shown in Fig. |2h, is equally satisfactory. The most interesting values from the fits are of 
course those for DceCu6, which are much lower than those for Cu. Assuming a Fermi velocity 
of the order of 10 3 m/s, a value for DceCue of 0.1 x 10~ 4 m 2 /s yields a mean free path l e 
of 3 nm, which is not surprising in view of the strongly granular nature of the films. Still, a 
low value for Dn by itself does not necessarily reflect the heavy Fermion character : if CeCug 
is taken as a Cu matrix strongly diluted by a small amount of Ce atoms, with a mean free 
path of the order of the interatomic distance, D would also be very small, of order 10~ 4 m/s. 
However, as we show now, just a low value for cannot describe the measurements. To 
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demonstrate this, we calculated T c (ds) for Nb/CeCu6 with the parameters given above, but 
with the DOS-ratio for Nb/Cu, thereby mimicking a very dirty but otherwise standard N- 
metal. The result, shown in Fig. [2h as dotted line, shows that T c now drops at much lower 
d,Nb- This can be understood by realizing that the low diffusion constant inhibits penetration 
of Cooper pairs in the N-metal and therefore also inhibits pair breaking, leading to a smaller 
amount of suppression of T c . What makes the difference is the high DOS-value for CeCu6 : 
the large number of available states works as a sink for Cooper pairs which counteracts the low 
diffusion constant. Our major conclusion is therefore that both the low Dn and the high Nn 
are necessary ingredients in the description of the data. This leads to the question what the 
effective mass needs to be in order to suppress the DCO which is so characteristic for an S /N 
system. For this we calculated values for dT/dH^Tc, the inverse parallel critical field slope 
at T c , as function of the effective mass m* /m e (m e being the bare electron mass) as used in 
the free-electron expressions for Nx = and Dx = , with kp the Fermi wave vector. 

As shown in Fig. low values for m* yield a finite value for the inverse slope, signifying 3D 
behavior and therefore a DCO, which goes to around m* = 10, meaning that the DCO has 
disappeared and 2D behavior has set in. Clearly, CeCu6 is far into this regime. 

Finally, we come back to the difference in proximity effects between F-, and HF-metals. 
Using the expression for £jy given in the Introduction at a typical value of T = 5 K and with 
the fitted values for DceCuQ, we find S,ceCu6 ~ 1-5 nm. This is very similar to values found 
for strong ferromagnets, but the physics of the strong suppression of superconductivity which 
is found both in the F- and HF-systems is different. In the F-case, the small value of £p 
derives from the large pair breaker (E ex ), in the HF-case from the small Dx- As shown in 
Fig. a low value for Dx does not yield strong suppression of superconductivity, that is 
actually due to the high value for Nx- This emphasizes once more that the basic proximity- 
effect parameter is not £f,hf but rather 7 = (ps£,s)/(px£,x) with X = F, HF (3"HTT|. Since 
px oc (NxDx) -1 and £x oc \JTJx it can easily be seen that for the F-system 7 is large due 
to E ex , and for the HF-system 7 is large due to NceCu6- For the HF-system, an extra pair 
breaker is not needed in the description. 

* * * 

This work is part of the research program of the 'Stichting voor Fundamenteel Onderzoek 
der Materie (FOM)', which is financially supported by NWO. We would like to thank R. 
Boogaard for preparing and measuring the Cu/Nb bilayers. A.O. acknowledges support from 
the ESF program 'FERLIN'. 



REFERENCES 
[1] de Gennes P. G. , Rev. Mod. Phys. 36 (1964) 225 . 

[2] Miihge Th. M., Garif'yanov N. N., Goryunov Yu. V., Khaliullin G. G., Tagirov L. R., Westerholt 

K., Garifullin I. A. and Zabel H., Phys. Rev. Lett. 77 (1996) 1857 . 
[3] Aarts J., Geers J. M. E., Briick E., Golubov A. A. and Coehoorn R., Phys. Rev. B 56 (1997) 

2779. 

[4] Hewson A. C, in 'The Kondo problem to Heavy Fermions', Cambridge Studies in magnetism, 

part 2 (Cambridge University Press, 1993). 
[5] Nieuwenhuys G. J., in 'Handbook of Magnetic Materials', ed. K. H. J. Buschow, Volume 9, pi 

(Elsevier Science, 1995). 



A. OTOP et al: PROXIMITY EFFECTS IN NB / CeCU 6 



7 



[6] Schlager H. G., Schroder A., Welsch M. and von Lohneisen H., J. Low Temp. Phys. 90 (1993) 
181 . 

[7] Banerjee I. and Schuler I. K., J. Low Temp. Phys. 54, (1984) 501. 

[8] Chun C. S. L., Zheng G.-G., Vincent J. L. and Schuller I. K., Phys. Rev. B 29 (1984) 4915. 
[9] Takahashi S. and Tachiki M., Phys. Rev. B 33 (1986) 4620. 
[10] Ciuhu C. and Lodder A., Phys. Rev. B 64 (2001) 224526. 

[11] Groten D., van Baarle G. J. C, Aarts J., Nieuwenhuys G. J., and Mydosh J. A., Phys. Rev. B 
64 (2001) 144425. 

[12] Amato A., Jaccard D. and Walker E., Sol. St. Comm. 58, 507 (1985). 

[13] Geers J. M. E., Hesselberth M. B. E., Aarts J. and Golubov A. A., Phys. Rev. B 64 (2001) 
094506. 

[14] Park S. I. and Geballe T. H., Physica 135B, 108 (1985). 

[15] Banerjee I., Yang Q. S., Falco C. M. and Schuller I. K., Solid St. Comm. 41 (1982) 805. 

[16] The values for Rb given in ref. |10| are quoted in units fiflcm and were obtained by somewhat 
arbitrarily multiplying the fitted value for Rb by a sheet thickness of 10 nm. A value such as 
3.17 fj,Qcm in ref. |10| therefore corresponds to 317 10 -8 fj,Qcm 2 in the present paper. 

[17] A. A. Golubov, Proc. SPIE 2157 (1994) 353. 



